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HORIZONTAL TAIL LOADS IN MANEUVERING FLIGHT 1 


By Hexby A. Peabsox, Wiluam A. McGowax, and James J. Donegan 


SUMMARY 

A method is given for determining the horizontal tail loads in 
maneuvering flight. The method is based upon the assignment 
of a load-factor variation with time and the determination of a 
minimum time to reach peak load factor. The tail load is 
separated into various components. Examination of these com- 
ponents indicated that one of the components was so small that 
it could be neglected for most conventional airplanes; therefore, 
the number of aerodynamic parameters needed in this computa- 
tion of tail loads was reduced to a minimum. 

In order to illustrate the method, as well as to show the effect 
of the main variables, a number of examples are given. 

Some discussion is given regarding the determination of maxi- 
mum tail loads, maximum pitching accelerations, and maximum 
pitching velocities obtainable. 

INTRODUCTION 

The subject of maneuvering tail loads has received con- 
siderable attention both experimentally and theoretically. 
Theoretically, methods and solutions have been derived for 
determining the horizontal tail load following either a pre- 
scribed elevator motion (references 1 to 3) or an assigned 
load-factor variation (reference 4). 

The first approach has been adopted into some of the load 
requirements where the type of elevator movement specified 
consists of linear segments whose magnitudes and rates of 
movement are governed by the assignment of a maximum 
initial elevator movement consistent with the pilot’s strength. 
The rates of movement and the time the elevator is held 
before reversing are so adjusted that the design load factor 
will not be exceeded. 

The results of reference 5 show, as is to be expected, that 
only when the aerodynamic force coefficients are accurately 
known from wind-tunnel tests can good agreement be ob- 
tained between measured and calculated tail loads. At the 
design stage, however, only general aerodynamic and geo- 
metric quantities are available and some of the more im- 
portant stability parameters are not known accurately. 
Thus, the work involved in the solution for the tail load 
following a given elevator motion is not considered to be in 
keeping with the accuracy of the results obtained. Con- 
sequently, there appears to be a need for an abbreviated 


design method of computing tail loads which, although 
incorporating approximations, will nevertheless be based on 
the theoretical considerations of the problem. 

If the load-factor variation with time is specified and the 
corresponding tail load, elevator angles, and load distribu- 
tions are subsequently determined, a simpler and equally 
rational approach to the tail-load problem can be made. 
Although this approach has been used, to a limited degree 
(reference 4), several shortcomings have limited its use. 

The purpose of this report is to develop further the load- 
factor or inverse approach and to present a method of com- 
puting horizontal tail loads which is comprehensive and 
generally simple. To this end, (1) the shape of the load- 
factor curve and the minimum time required to reach the 
peak load factor have been determined from an analysis of 
pull-up maneuvers that were available, (2) the minimum 
time required to reach the peak load factor has been deter- 
mined from a theoretical analysis which is supported in some 
measure by statistical data obtained from a number of 
flight tests with airplanes of widely varying sizes, and (3) 
the equations relating the various quantities are presented. 
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SYMBOLS 

wing span, feet; shape factor in equation (13) 
tail span, feet 
chord, feet 

mean aerodynamic wing chord, feet 
lift coefficient ( L/qS) 

pitching-moment coefficient of airplane with- 
out horizontal tail ( Mb/gS *) 
pitching-moment coefficient of isolated 
horizontal-tail surface 

acceleration due to gravity, feet per second 
per second 

pitching moment of inertia, slug-feet* 
radius of gyration about pitching axis, feet 
empirical constant denoting ratio of damping 
moment of complete airplane to damping 
moment of tail alone 
Lift, pounds 

local lift at any spanwise station 
airplane mass, slugs (W/g) 
pitching moment, foot-pounds 
airplane load factor at any instant 


1 Supersedes NACA TN 2U7S, “Horizontal Tall Loads in Maneuvering Flight” by Henry A. rearson, William A. McGowan, and James J. Donegan, 19S0. 
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maximum increment in load factor 
dynamic pressure, pounds per square foot 


(y v ‘) 


5 wing area, square feet 

S t horizontal-tail area, square feet 

t time, seconds 

U time to reach peak of elevator deflection, 

seconds 

V airplane true velocity, feet per second 

W airplane weight, pounds 

x t length from center of gravity of airplane to 

aerodynamic center of tail (positive for con- 
ventional airplanes), feet 
y spanwise dimension, feet 

y* nondimensional spanwise dimension 

a, b, e )constants occurring in equations (13), (23), 

A, B, O, D, E) (26), and (30) 

Kj, K tl Kt constants occurring in basic differential equa- 
tion (see equation (3)) 

X time to reach peak load factor, seconds 

p mass density of air, slugs per cubic foot 

17 ( tail efficiency factor (g_t/g) 

a wing angle of attack, radians 

a average angle of attack of horizontal stabi- 

lizer, radians 

a, tail angle of attack, radians 

/3 angle of sideslip, degrees 

7 flight-path angle, radians 

6 attitude angle, radians (a+7) 

5 elevator angle, radians 

e downwash angle, radians aj 

i t tail setting, radians 


The notations a and 6, a and 6, and so forth, denote 
single and double differentiations with respect to t. 

The symbol A represents an increment from the steady- 
flight datum value. 

Subscripts: 

0 initial or selected value 

t tail 

max maximum value 

l o zero lift 

geo geometric 

c camber 

METHODS 


METHOD OF DETERMINING THE DYNAMIC TAEL LOAD 

Basic equations of motion. — The simple differential equa- 
tions for the longitudinal motion of an airplane for any 


elevator deflection (see method given in reference 2) may 
be written as 

A “ qS ~(irf) t Vt * St A5=0 (1) 

2 ^ — AL t x t -\ — AS — mk y i 8=0 (2) 

Equations (1) and (2) represent summations of forces 
perpendicular to the relative wind and of moments about the 
center of gravity. (See fig. 1 for direction of positive quanti- 
ties.) Implicit in these equations are the following assump- 
tions: 

(1) In the interval between the start of the maneuver and 
the attainment of maximum loads, the flight-path angle does 
not change materially; therefore, the change in load factor 
due to flight-path change is small. 

(2) At the Mach number for which computations are 
made, the aerodynamic derivatives are linear with angle of 
attack and elevator angle. 

(3) The variation of speed during the maneuver may be 
neglected. 

(4) Unsteady lift effects may be neglected. 

By use of the relations 8=y+a, 8=y+a, and 0=7+a, 
equations (1) and (2) are reducible to the equivalent second- 
order differential equation 

a-\-K\a-\-Kt A a — Kz AS (3) 

where 


pVVdC^StX? (K .rf.VrfCac"! 

Kt -*z L^r ts " \^ t + d- a r~d ^ 6 J 

T , pV* [dC m S’ dC Lt s t x,rs, de\, 

2m { da ky'b da t Vt IV L.V d a)* 

dC L K^ L SxT\) 

^2 m J| 

and 

K _ P v y s.*. , 

Ki ~2m\ dd Vt + 

dC m , S t * dC Ll dC Lt KV p VSA 

dS 1,1 b t k Y * da t dS 2 mk Y 2 ) 

In equations (1) and (2), Aa, 7, 6, AS, and A L, will, in a 
given maneuver, vary with time. Using the relations between 
6, 7, «, and their derivatives permits equation (2) to be re- 
written as follows to give the increment in tail load: 

„ T dC n * S* mk Y *a mk Y *7 , dC m , S t * Ar fAS 
h. — x, (4) 
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In a still shorter form, equation (4) may be written as 

AL [ =ALt a J rAL t;i J t-AL t --{-AL to (5) 

Equations (4) and (5) show that the tail-load increment (the 
increment above the steady-flight datum value) at any time 
is composed of four parts: AL ta , associated with the angle-of- 
attack change; AL t _, associated with angular acceleration 
about the flight path; A L,^, associated with angular accelera- 
tion of the flight path; and AL t<; , required to compensate for 
themoment introduced by change in camber of the horizontal- 
tail surface. The load A L tc is generally small but in some 
extreme configurations may amount to 10 percent of the 
total increment and thus, for the present, it is retained in 
the development. 

If the load-factor-increment variation with time An is 
known, then by the usual definition 

<«> 


and 


a 


hW/S 
dC L 


da 


2 


The following relation also exists between An and 7 : 

A ng=yV (7) 


so that 


z=M 
7 V 


( 8 ) 


When equations ( 6 ) and ( 8 ) are substituted into equations 
(4) and (5), the four tail-load components then become 


. . iO. It’S 

(9a) 

***- „ iO L n 

gS^ da 

(9b) 

Wkr 2 . 

(9c) 

Ar dC mt S? AS 

d5 Vtg - b t x t 

(9d) 


Thus, if the variation of the load factor with time An and 
the geometric and aerodynamic characteristics of the airplane 
were known, the first three components of the tail load could 
be found immediately. The magnitude of the fourth com- 
ponent, that due to horizontal-tail camber, would follow 
from equation (3) in which the elevator angle is seen to be 


. . a . Ki . - K 2 * 
A5 =j=-+^ra+jf - A a 

J-A.3 


( 10 ) 


Substitution into equation (10) of the values of Aa, a, and 
a from equation ( 6 ) yields the value of the elevator angle at 
any instant 


A S = 


W/S 




An) 


(ID 


so that, finally, the fourth component is given as 


so that 


Aa c 


An W/S 

da 2 


AW/S 

z dC L 

da ® 


da 


SfWJS,- 




da 


Ki 


The procedure outlined shows that the tail-load magnitude 
can be determined if the load-factor variation is known. 

Types of load-factor variation. — The relation between the 
tail load, the geometric and aerodynamic characteristics, 
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and the load factor having been established, it is desirable to 
establish a load-factor variation which is reasonable as well 
as critical insofar as loads are concerned. The maximum 
value of load factor is usually specified; however, there are 
many possible variations for the shape. Regardless of the 
details of shape, the load factor may be considered to rise 
smoothly and continuously to a maximum, the rate of rise 
depending upon several variables. Beyond the maximum 
value of the load factor the return to initial conditions can, 
at the will of the pilot, be either gradual or rapid. 

Experiments as well as theoretical studies have already 
indicated that the maneuver that combines maximum 
angular and linear accelerations causes critical loads in both 
the wing and tail. One such maneuver occurs when the 
maximum load factor is reached as rapidly as possible by 
using an initial elevator movement which is greater than 
that required to reach a given steady-trim value of the load 
factor. This initial elevator movement is followed by a 
rapid checking of the maneuver either by returning the 
elevator quickly to neutral or by reversing the controls. 

The shape of the load-factor curve for such a maneuver 
may be expressed approximately by several analytic func- 
tions, one of which is 

A n^afe- 01 (13) 

By way of illustration, figure 2 shows details of the shape of 
the load-factor curve obtained with the use of equation (13) 
for which the constants have been adjusted so that an 8g 
peak is reached in 1 second. By further adjustment of the 
constants the load factor can, within certain limits, be made 
to rise to any specified peak and to diminish in any prescribed 
manner. 

Because the positive slopes obtained from equation (13) 
are always greater than the negative slopes, the positive 
angular accelerations are greater than the negative ones. In 
general, this condition is true for most high g critical ma- 
neuvers performed by most classes of airplanes, but maneuvers 
may occasionally be performed for which the reverse may 
be true, particularly for small airplanes. 



Determination, of constants. — From equations (9), (11), 
and (12) the required quantities relating to load factor are 
seen to be An, h, and fi. Since the increment An is to be 
given by 

An=at*e~ ct (13) 

then at maximum load factor 

n=0=An(j—c' S J (14) 

Thus t=— at maximum load factor. Let N 

(!)“'- 

so that 

AW_L\V-“ 

N~\bfc) e 

Let — =X. Then 

#-G) V( " D 

Equation (17) is in nondimensional form where X is the time 
to reach the peak load factor and b is a constant. 

When equation (17) is differentiated, the first and second 
derivatives become 


fix An , ( 1 . \ 

W = N 6 U/X~V 

00 

V — / 

n\’ 4»„rx’/, 1 \ 2 X ..1 

(19) 


In equations (17) to (19) the quantities A r , X, and b are 
now required in order to determine the variation of An, fi, 
and fi. The value of N is immediately available from the 
required maneuver load factor; whereas the time to roach 
the peak load factor X can be obtained from examination of 
available records or by specification. The constant b, as 
may be seen from equation (17), can best be described as a 
“shape” factor and has no particular physical significance. 

The values of X and b should be associated with a maneuver 
which produces maximum tail loads; therefore, the time 
X to reach peak load factor should be the minimum possible 
consistent. with possible pilot action and airplane response. 
The shape factor b should also be consistent with both of 
these. 

In connection with the determination of the minimum time 
to reach peak load factor, the results shown in figure 3 for 
a typical airplane are informative. Figure 3 (a) shows the 
load-factor variation following several abrupt jump elevator 
movements. The load factor varies with the elevator posi- 
tion, but the time to reach peak load factor does not. Figure 
3 (b) shows the load-factor variation for several abrupt hat- 
shape elevator impulses. Again the load factor is seen to 


= A n max . Then 

(15) 

(16) 
(17) 
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(a; Jump. (b) Impulse. 

FlGUKi 3.— Incremental-toad factor variations following control movement. 



vary with the amount of elevator deflection but the time to 
reach the peak value remains constant. Although the time 
to reach the peak load factor shown in figure 3 (b) remains 
constant, it is seen to be less than that shown in the previous 
case; therefore, an impulse elevator motion produces a smaller 
value of X than the jump type. 

Because of inertia and elasticity in the control system, the 
pilot cannot move the elevator instantaneously but requires 
some finite time f t to do so. A possible critical type of eleva- 
tor impulse thus appears to be one which increases linearly 
to maximum and decreases at the same rate to zero. In order 
to determine the minimum time to reach peak load factor 
associated with such a variation, the equation of motion 
(equation (3)) has been solved for the tri ang ular elevator 
impulse for airplanes of various static stabilities and damping. 

The results of the computations are given in figure 4 in 
which the minimum time X to reach peak load factor is 
plotted against the time fi required to deflect the elevator. 

For completeness the curves of figure 4 are labeled for the 
actual values of K a employed in the computation as well 
as for relative values of stability. By a series of computa- 
tions the damping term Ki was found, as was to be expected, 

213eST— 63 9 


to have only a secondary effect on X. The curves apply to 
an average value of the damping constant. The upper 
curve, labeled “low stability,” should be associated with 
rearward center-of-gravity positions (that is, low static 
margin) in combination with one or both of the following: 
low dynamic pressure or heavy airplanes. The lower curve, 
labeled “high stability,” would be associated with forward 
center-of-gravity positions in combination with one or both 
of the following: high dynamic pressure or light airplanes. 
It is seen that X increases almost linearly with t L and also 
increases when the restoring forces are reduced, that is, 
when the stability is reduced. 

A preliminary value of the shape factor b (required in 
equations (17) to (19)) was initially determined from flight 
records of typical impulse maneuvers by fitting curves of 
the type given by equation (13) through several points of 
the actual time histories and determining the constants. 
The results of this first step were then modified, by the 
results of the same computations which had been made to 
determine X, and the variation of b with fi given in figure 5 
was obtained. Since the b factor is not found to be critical, 
an average value of 5.0 is suggested, although as a refinement 
the values from figure 5 may be used. 
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Fiocre 5.— Variation of shape factor 6 with h. 

The question of the value of U to use is one which must be 
solved either from experience or from a knowledge of the 
characteristics of the controls and the control system. For 
conventional airplanes having the usual amounts of boost 
and no rate restrictors, the following values of ti are sug- 
gested as representative: 

h 

Fighters or small civil airplanes with weight limit from about 500 

to 12,000 pounds, seconds 0. 20 

Two-engine airplanes with weight limit from 25,000 to 45,000 

pounds, seoonds 0. 25 . 

Four-engine airplanes with weight limit from 50,000 to 80,000 

pounds, seconds 0. 30 

Airplanes with weight limit above 100,000 pounds, seconds 0. 40 

The mini mum time X given in figure 4 was actually estab- 
lished separately from the adopted load-factor variation; 
therefore, in applying the inverse method, the derived ele- 
vator impulse would not be expected to agree in detail with 
the “tent” type impulse used in the derivation. 

The first three tail-load components can now be computed 
by inserting the values of An, n, and n from equations (17) 
to (19) into equation (9) and using appropriate values of X 
from figure 4. In order to facilitate this computation, 
curves of An /N, nk/N, and nk 2 /N plotted against tj\ are 
given in figure 6 for the suggested value of 6=5. Actually, 
in order to apply the results of figure 6 it is convenient to 
find first the components A L ta , and so forth, in terms of the 
nondimensional time f/X and then to convert to time t in 
seconds. In order either to compute the fourth component 
or to obtain the elevator angles for use in chord loading, the 
constants K u 1 Kj, and K a of equation (3) must also be known. 

Thus, in terms of f/X and the ordinates of figure 6, the 
various tail-load components are 

AL ‘ a ~rfr~ (° rdinate of %• 6 ( a )) (20a) 


^ (Ordinate of fig. 0(c)) (20b) 

ci OtC l A 

JSJ*. -J- 

Wl. 2 AT 

A L tx . = — y^~ y (Ordinate of fig. 6(b)) (20c) 

A r d.C m , S t * W N ["Ordinate of fig. 6(c) , 

tc = ~dT Vt Wt ~s^7dC, L 3? + 

K> da 

K, (Ordinate of %. W> +g|(0rdimlg of ^ 6(a)) j 

(20d) 

The constants K u K % , and K 3 defined previously herein are 
the same as those given in reference 2, except for changed 
signs caused by specifying x, as positive. 

The conversion to time f is made by multiplying values of 
the base scale f/X by X. 

Sample calculations for incremental tail loads. — The 
results of several examples are given to illustrate not only the 
method but also the effect of each of a number of variables 
on the incremental tail load of a fighter-type airplane, the 
geometric and aerodynamic characteristics of which are 
given in table I. In order to illustrate the effect of static 
stability, results have been computed for three center-of- 
gravity positions with the assumption that an 8 g recovery is 
made at 19,100 feet from a vertical dive at an airspeed of 400 
miles per hour. In order to illustrate the effect of the time 
of the elevator impulse on the tail load, computations were 
carried out at one of the center-of-gravity positions for 
several values of fj. The cases considered and the airplane 
characteristics are given in table I. 

TABLE I.— AIRPLANE CHARACTERISTICS 
(a) Geometric. 


Gross wing area, S, square feet 300 

Gross horizontal-tail area, S t , square feet 60 

Airplane weight, W, pounds 12, 000 

Wing span, 6, feet _ 41 

Tail span, 6 ( , feet 10 

Radius of gyration, hr, feet 6. 4 

Distance from aerodynamic center of airplane less tail to 
aerodynamic center of tail, x t , feet: 

Center of gravity, 29 percent M. A. C 20. 0 

Center of gravity, 24 percent M. A. C 20. 3 

Center. of gravity at aerodynamic center.. 21. 0 

(b) Aerodynamic. 

Slope of airplane lift curve, dCJda, radians , 4. 87 

Slope of tail lift curve, dChJdat, radians 3. 15 

Downwash factor, dtjda 0. 54 

Tail efficiency factor (?</?), vt X 00 

Empirical airplane damping factor, K L 1 

Elevator effectiveness factor, dCiJdS, radians 1. 89 

Rate of change of tail moment with camber due to elevator 

angle, dC m Jd5, radian —0. 57 

Rate of change of moment coefficient with angle of attack 
for airplane less tail, dC m /da, radians: 

Center of gravity, 29 percent M. A. C. 0. 625 

Center of gravity, 24 percent M. A. C 0. 403 

Center of gravity at aerodynamic center 0. 000 




AA/N 
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Florae 8.— Variation of fneremental-load-&ctor curves wltli time ratio. 
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The specified conditions for the sample computations are 
given in table II. The computed results for tail components 
are given in figures 7 and 8. Figure 7 gives results for vary- 
ing the center of gravity and figure 8 gives similar results for 
varying ti. The tail-load components are computed from 
equations (20) and the derived elevator angles from equation 


(11). If the increment in tail load due to camber and the 
incremental elevator angle are not required, the values of K 
need not be computed and the computations are considerably 
shortened. Figures 7 and 8 show that a maximum error of 
only about 4 percent is introduced by this omission. 



Piques 7.— Effect of center -of-gravity position on Incremental-tall-load components. 
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TABLE II.— SPECIFIED CONDITIONS OF SAMPLE 


PROBLEM 

Increment in load factor 8. 0 

Altitude, feet 19, 100 

Air density, slug per cubic foot 0.001306 


Case 

c- ff- 

(percent 

Ml.C.) 

ft 

Kx 

Hi 

Ex 

(flg.ti 

1 

a. c. 

0.2 

4.93 

30.4 

-33.4 

0.45 

2 

24 

.2 

4.72 

15.2 

-32.2 

.50 

3 

2D 

.2 

4.61 

8.45 

— 3L 7 

.56 

4 

24 

.4 

4.72 

15.2 

-32.2 

. 1 1 

5 

21 

.6 

4.72 

16.2 

-32.2 

L 02 


METHOD OF DETERMINING THE TOTAL TAIL LOAD 

The initial or steady-flight tail load and elevator angles to 
which the computed incremental values are to be added 


must also be determined. In steady flight, the horizontal 
tail furnishes the moment required to balance the moments 
from all other parts of the airplane so that the initial load 
may be written as 




C mo zSc 

x t 


dC m WS 
dG L bx t 


cos Yo 


( 21 ) 


Thus the total tail load at any time in a maneuver is com- 
posed of the four previously mentioned parts plus the 
components given in equation (21). Only the first term of 
equation (21) represents a new type of load because the 
second term is a load of the type given by equation (9a) or 
equation (20a) and its effect may be immediately included in 
the computations by multiplying the ordinates of figure 6 by 
N + cos to instead of by A T . 




77/7!e, sec 

Fig cm 8.— Effect of ft on fmzemental-tatl-load components. 


Time, sec 
Center o f gravity, 0.24J. 
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The initial elevator angle required to balance the airplane 
in steady flight varies with airplane C L and center-of- 
gravity position so that, in general, 5 0 must be obtained from 
wind-tunnel data. Without results of wind-tunnel tests, a 
rough rule which can be used as a guide at the design stage 
in determining the elevator position is that the final elevator 
setting will be so adjusted by repositioning of the stabilizer 
setting during acceptance tests that it will be near a zero 
position at the cruising speed and at the most prevalent 
center-of-gravity position. 


METHOD OF DETERMINING MAXIMUM VALUES 

Maximum tail loads and angular accelerations. — The 
method outlined enables a point-by-point evaluation to be 
made of the quantities that determine the tail load. Such 
detail may often be unnecessary and the procedure may be 
shortened by evaluating only those points near the load 
peaks or, alternatively, by accepting an approximation to 
the results. One such approximation which may be made 
is to balance the airplane at the combinations of load factor 
and angular acceleration which would result in maximum 
up and down tail loads. 

Figure 7 shows that the maximum down tail load in a 
pull-up occurs near the start of the maneuver and before 
appreciable load factor is reached. This maximum load is 
practically coincident with the negative maximum in the 
L t s tail-load component. Since, for a given configuration, 
this component increases as the center of gravity is moved 
forward and since the steady-flight down load increases with 
speed, the maximum down tail load in a pull-up occurs at 
the highest design speed in combination with the most 
forward center-of-gravity position. 

Figures 7 and 8 show that at the time of the maximum 
down-tail-load increment the elevator is near but has not 
quite reached its peak position. Also at the time of maxi- 
mum up-tail-load increment the elevator is near its zero 
position, although it may be on either side of this position 
depending upon the stability and the time These results 
suggest that the maximum down load for the elevator and 
the hinge brackets would occur with the airplane center of 
gravity well forward and at the start of the maneuver. 
The maximum load for the stabilizer is likely to occur at 
the peak load factor. 

Figure 7 also shows that the up tail load occurs near the 
peak of the component L ta as well as near the positive maxi- 
mum peak in the component L t ~. Since the component 
L, a increases as the center of gravity is moved rearward 
and since a decrease in speed generally reduces the initial 
down load, the maximum up tail load occurs at the upper 
left-hand corner of the V-n diagram for the most rearward 
center-of-gravity position. 

The maximum tail load in a pull-up maneuver may be 
written as 


Lt. 


Cm 0 jSc dC m 

x, ~ t " d,C L 


WS , , .16 

( 22 ) 


the load-factor increment at maximum down load is nearly 
zero and at maximum up load it is nearly equal to N so 
that if the positive and negative values of 0 max can bo de- 
termined, a relatively simple method for determining maxi- 
mum loads is available. 

Since by definition e=a+y, an expression for angular 
acceleration can be derived from equations (6) and (7) and 
written in the form 


by 


s W/S-,g . 

The maximum angular acceleration can be approximated 


max dC L V\ 

da ® 


(23) 


For the maximum positive pitching acceleration, B is the 
maximum positive ordinate in figure 6 (c) and C is the ordi- 
nate of figure 6 (b) at a value of tj A for which B was deter- 
mined. Thus, B is 6.5 and C is 0.95 for this example. 

For the maximum negative pitching acceleration, B is the 
maximum negative ordinate in figure 6(c) and O is the 
ordinate of figure 6(b) at a value of t/\ for which B was 
determined. Thus, B is — 5.8 and C is 0.80. For use in 
equation (23) the values of X for the maneuver are available 
from figure 4 and the other quantities are available from 
the conditions of the problem. The maximum loads can 
be given by the following equations: 

For maximum up tail load in the pull-up: 


r.. C mo qSc , dC m TVS' 

‘ ntax+ “r* * ~ t 1 


x, 


dC L bx, 




(24a) 


For maximum down tail load in the pull-up: 




x t dCi, bx, 


x, 


XI dCt 

\ da 


q\ 


V 
(24b) 


) 


For push-downs to limit load factor, equations (24a) and 
(24b) still apply with changed signs for N and changed 
directions for and Z w _. A question arises as to 

whether the maximum down tail load at the start of a pull-up 
with forward center-of-gravity position is greater than that 
which would occur when pulling up from a negative load- 
factor condition with the center of gravity in the most 
rearward position. This can be determined only by com- 
puting both cases and seeing which is the larger. 

Maximum value of angular velocity. — The maximum 
value of the pitching angular velocity in the pull-up may 
also be found in a manner similar to that used to obtain the 
maximum angular acceleration. Since 6=ct-j-y and the 
relations involving these quantities in terms of load factor 


where the sum of the second and third terms is to be a 
maximum in the maneuver. From the previous discussion 
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are given by equations (6) and (7), the following equation 
may be written: 

* nW/S,Ang 

6 -dor + ~v~ f25) 

da ® 

The maximum angular velocity may be approximated by 




where D is the maximum positive ordinate in figure 6 (b) 
and E is the ordinate of figure 6 (a) at a value of f/X for 
which D was determined. Thus D, for this example, is 1.95 
and E is 0.48. 

In the steady turn or pull-up at constant g, the angular 

q N 

velocity is usually given by the expression $=1.0 The 

difference between the factor 1.0 of this expression and the 
factor 0.48 of equation (26) is more than made up by the 
angle-of-attack component of the angular velocity. 

APPBO XIMATE METHOD OF DETERMINING LOAD DISTBIBUTION 

Symmetrical loading. — The spanwise distribution of the 
total load can be formulated with various degrees of exact- 
ness. If information regarding details of the angle-of-attack 
distribution across the span were available, then an exact 
solution could be obtained for the loading with the use of 
existing lifting-surface methods. The following method may 
be used as a first approximation to the solution. 

From the total tail load, the total tail lift coefficient C L , 
can readily be found. The average effective angle of attack 
a of the stabilizer portion is given in the definition 

Cl, = c [a 5 1 dy*+ c,,(5 0 +A5) | dy * (27) 

where only a is assumed as unknown and c ta and c if may be 
taken as the rates of change of section lift coefficient with a 
and 5, respectively. 

Thus, for constant elevator angle across the span, 


Cl, — (5o+A.5) J* Cjj = dy * 

S " J> 4^ 


In a practical case both integrals in equation (28) need be 
evaluated only once for a given configuration and Mach 
number. A plot of a against C L , with 5 as a parameter 
would be useful in further computations. With a known as 
a function of Cl c and 5, the local lift at any spanwise station 
is then obtained from the expression 

l=c h qc=[c la a+c lt (8o+M)]gc (29) 


TJnsymmetrical loading. — Up to this point the total loads 
have been assumed to be symmetrical about the airplane 
center line, whereas, in reality, the load may have an un- 
symmetncal part. The sources of this dissymmetry may 
be due to uneven rigging, differences in elasticity between 
the two sides, or to effects of slipstream, rolling, and sideslip. 
The first two sources are usually inadvertent ones while the 
last two are difficult to determine without either wind-tunnel 
tests or a knowledge of how the airplane will be operated. 
Present design rules regarding dissymmetry of tail load are 
concerned more with providing adequate design conditions 
for the rear of the fuselage than with recognizing that at the 
m ax im um critical tail load some dissymmetry may exist. 

Tests in the Langley full-scale tunnel (reference 6) and 
flight tests (reference 7) of a fighter-type airplane, as well 
as unpublished flight tests of another fighter-type airplane, 
indicate that the tad-load dissymmetry varies linearly with 
angle of sideslip so that the difference in lift coefficient 
between the two sides of the tail can be given as 

Cl. -C Lt (30) 

‘Right ‘Lett 

The average values of A per degree found for the two fighter- 
type airplanes are approximately 0.01. No similar values 
are available for larger airplanes nor for tail surfaces having 
appreciable dihedral. 

In maneuvers of the type considered herein it is doubtful 
that angles of sideslip larger than 3° would be developed at 
the time the maximum tail load is reached. If the value 
of the sideslip angle at the time of maximum tail load can 
be established, equations (27) to (29) are easily modified to 
include this effect, provided that the approximate value of 
A is known. 

Chordwise loading. — The chordwise distribution can be 
determined for any one spanwise station in either of two ways. 
One way for design work is outlined in reference 8. A 
knowledge of the airfoil section and the quantities contained 
in equation (29) suffices for this determination. 

If pressure-distribution data are available for a similar 
section with flaps, an alternative way would be to distribute 
the load chordwise according to the two-dimensional pressure 
diagrams with the use of the computed values of section lift 
coefficient and elevator angle. 

DISCUSSION 

The method presented is another approach to the deter- 
mination of tail loads. From the results given in figures 7 
and 8, it can be seen that the camber component L te is so 
small that for all practical cases it may be omitted with 
considerable simplification in the computation of tail loads. 
This omission reduces to a minimum the number of aero- 
dynamic parameters needed to compute the tail loads. 
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It is possible, in the application of the present method with 
the use of the suggested values of t u that the derived elevator 
angles may not be within the pilot’s capabilities. Since it 
must be assumed that all airplanes, to be satisfactory, 
should have sufficient control to reach their design load 
boundaries, such an occurrence requires only that the time 
to reach elevator peak deflection h be increased so as 
to reduce the elevator angle. The results of figure 8, in 
which the time 4 is varied, furnish a usoful guide for deter- 
mining the increase U that might be required. 

If sufficient information is available, it is recommended 
that existing lifting-surface methods be used in determining 
the spanwise distribution of the total load; however, if 
information of the angle-of-attack distribution across the 
span is not known, the method presented may be used as 
a first approximation. 

Along some of the boundaries of the V-n diagram tail 
buffeting may occur. Measurements show that buffeting 
usually occurs along the line of ma xim um lift coefficient and 
again along a high-speed buffet line which is associated with 
a compressibility or force break on some major part of the 
airplane. All airplanes are subject to buffeting at the 
design conditions associated with the left-hand comer of the 
V-n diagram. Only high-speed or high-altitude airplanes 
or both are capable of reaching the other boundary. Meas- 
urements show that the oscillatory buffeting loads may be 
so high that the designer should at least be cognizant of 
them at the design stage. 

The maximum angular acceleration varies inversely with 
airspeed and directly with the load factor with the contri- 
bution due to acceleration in angle of attack likely to be 
more important than the angular acceleration of the flight 
path. A somewhat similar variation is indicated for the 
maximum angular velocity (equation (26)) where it is seen 
by direct substitution that the part due to angle of attack 
is likely to be larger than the part due to the angular velocity 
of the flight path. 

CONCLUDING REMARKS 

A simple method has been presented for determining the 
horizontal tail loads in maneuvering flight with the use of a 
prescribed incremental load-factor variation. 


The incremental tail load was separated into four com- 
ponents representing a, a, y, and c. The camber com- 
ponent L, e is so small that for most conventional airplanes 
it may be neglected; therefore, the number of aerodynamic 
parameters needed in this computation of tail loads was 
reduced to a minimum. 

An approximate method is presented for predicting maxi- 
mum angular accelerations and maximum angular velocities. 

The method indicates that maximum tail loads in a 
pull-up occur at forward centor-of-gravily positions and 
early in the maneuver. The maximum down tail loads in 
a pull-up occur at the highest design speed in combination 
with the most forward center-of-gravity position. The 
maximum up tail load occurs at the upper left-hand corner 
of the T r -n diagram for the most rearward center-of-gravity 
positions. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 9, 1950. 
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